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ABSTRACT 

The design and performance of a wide bandwidth hnear polarization modulator based on the Faraday 
effect is described. Faraday Rotation Modulators (FRMs) are solid-state polarization switches that 
are capable of modulation up to 10 kHz. Six FRMs were utilized during the 2006 observing season 
in the Background Imaging of Cosmic Extragalactic Polarization (BICEP) experiment; three FRMs 
were used at each of BICEP's 100 and 150 GHz frequency bands. The technology was verified through 
high signal-to-noise detection of Galactic polarization using two of the six FRMs during four observing 
runs in 2006. The features exhibit strong agreement with BICEP's measurements of the Galaxy using 
non-FRM pixels and with the Galactic polarization models. This marks the first detection of high 
signal-to- noise mm- wave celestial polarization using fast, active optical modulation. The performance 
of the FRMs during periods when they were not modulated was also analyzed and compared to results 
from BICEP's 43 pixels without FRMs. 



1. INTRODUCTION 

Numerous experiments in the last two decades 
have shown that the cosmic microwave background 
(CMB) is a powerful cosmological probe. The tem- 
perature anisotropy of the CMB has now been 
mapped to ex quisite prec is ion by many experi- 

ments (e.g; - Jones et al. I (120061^ : iReichardt et al. 

2009|) : iFried man e t al. I (I2OO90 : ,Sievers et al 



2009li: IJarosik et al. n[20Toh : IHlozek et al. I (1201^ 



IReichardt et al. I ()2011|) ). yielding tight constraints on 
the s ta ndard ACDM c o smolo gic al model (iBrown et al. I 
[20091: iLarson et al~l [2OIOI: iDunklev et al. I 120101: 
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IKeisler et an[20TlHSto"rv et al. lIMl) . 

Many current and upcoming experiments focus on 
measuring CMB polarization, which potentially encodes 
information from long before the epoch of matter- 
radiation decoupling. CMB polarization can be de- 
composed into two orthogonal bases: the "i?-mode", a 
curl-free mode, and the "B-mode", which possesses a 
curl or handedness. The shape of the E-mode power 
spectrum has been measured and is consistent with the 
ACDM cosmological model predicted by the temperature 
power spectrum. Multiple experiments have now mea- 



spectrum ( 


Leitch et al. 1 120051 iMontrov et al. 1 


2006; 


iPage et al. 


2007; Sievers et al. 112007: Wu et al. 1 


2007; 


iBischoff et al. 1 2008: Noha et al. 1 120081: Prvke et al. 



2009HChiang' et al. 1 120101: lOUIET Collaboration II2012D 
Primordial B-mode polarization, however, is not a 
facet of the standard ACDM model of the universe. In- 
stead, the inflationary paradigm predicts the existence 
of a stochastic gravitational wave background (GWB) 
that would imprint a unique pattern on the surface 
of last scattering (Polnarev 1985). While this faint 
signal does affect the CMB i?-mode and temperature 
power spectra, the GWB's effect is degenerate with 
other parameters. Only the i?-mode polarization pro- 
vides a unique and direct probe of the GWB. This sig- 
nal is expected to peak at degree angular scales with 
an amplitude dete rmined by the en ergy scale of infla- 
tion ( Ka mionkowskTand Kosowskv I [1997 ). A detection 
of the i?-mode signal would provide strong evidence for 
an inflationary phase in the early Universe. 

Detection of the sought-after B-mode signal presents 
many difficulties. First, the inflationary _B-mode ampli- 
tude is extremely small compared to sky temperature, 
the CMB's temperature, and even its temperature and 
£'-mode anisotropy. Currently, the most restrictive lim- 
its, derived from temperature anisotropy, put an upper 
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limi t on the expected i3-inode signal of less than 200 
nK (|Keisler et al. 11201 ItlStorv et al. II2012D . Second, po- 
larized emission from galactic and extra-galactic objects 
acts as a foreground contaminant for CMB polarization 
measurements. Even in the case of optimistic predic- 
tions for the signal amplitude, models suggest that fore- 
ground contamination could be an order of magnitude 
larger than the -B-mode sig nal at the largest angular 
scales (jAmblard et al. II20Q7D . Finally, exquisite control 
of systematic and instrumental effects, down to the tens 
of nK level, will be required before a detection of B- 
modes can be claimed with confidence. 

One approach to mitigating systematic errors is to 
modulate incoming polarization, thereby shifting the sig- 
nal to higher frequencies and away from lower frequency 
systematic contaminants. Rapid modulation eliminates 
concerns about time varying thermal, optical, or electri- 
cal fluctuations that often change on much longer time 
scales. In essence, modulation speed can be used to tune 
the signal band of the instrument, otherwise set by a 
combination of scan speed and beam size, placing the 
signal away from microphonic lines and low frequency 
1// noise. The lack of limitation on scan speed expands 
the parameter space for observation strategies. Modula- 
tion also mitigates polarized systematic effects that are 
introduced by optical elements between the modulator 
and the detector. 

There are additional benefits to polarization modu- 
lation. The reconstruction of the Stokes Q (or U) pa- 
rameter requires two independent measurements at dif- 
ferent detector angles, typically done by two detectors. 
Mismatched detector properties, such as those caused 
by caused by differential spectral response, differential 
pointing, and mismatched transfer functions, can result 
in spurious polarization. Modulation allows a single lin- 
early polarized detector to act as an independent po- 
larimeter, removing systematics resulting from combin- 
ing mismatched detectors. 

Many CMB polarimeters already employ mecha- 
nisms to modulate the incident radiation field about 
the optical axis of the instrument. The modulator 
effectively exchanges electric fields between two de- 
tectors (or a single detector in two orientations) and 
decomposes the radiation into two orthogonal bases. 
Traditional methods for modulation in millimeter wave 
polarimetry include physical rotatio n of the entire 
instrument abo ut th e optical axis (iHalverson et al. 
1998; Keating ] [200lt [P adin et al. 200 21: lYooiT'eTaL 
200& .O'SuUiyan et al. i .2008; .Hinderks et al. I 120091: 
Qgburnet al. I 120101) . rotation or translation of a 
wiregrid polarizer ([Caderni et al. I Il978t 'Chuss et al. I 
2012 ) or birefringent h al f-way ep late ([ Murray et al 
1992 IMdebrand et al. I [200(11: lOxlev et al. I l200 



Johnson et al. I 12007 
20101: iFilippini et al. 



lEssinger-Hileman et al. et al.l 

1201 It IKea~ et al. I 120111: 

Staniszewski et al. II2012I ). or parallactic angle rotation 
("sky rotation") with respect to the fixed instrument 
coordinate system. 

These classical modes of polarization modulation of- 
ten employ bulky and complex mechanisms, where any 
mechanical failure would result in the complete loss of 
polarization modulation. Faraday Rotation Modulators 
(FRMs) are coupled to individual detectors, eliminat- 
ing the single point failure mode and relieving the need 



for producing large grids and waveplates that are often 
limited to tens of centimeters in diameter. FRMs can 
be tuned individually for specific bandwidths, whereas 
birefringent crystal waveplates require specialized anti- 
reflection coatings that are difficult to optimize for mul- 
tiple frequencies. 

In addition, polarization modulation via mechanical 
motion is limited to a low frequency range due to toler- 
ances on vibrational and microphonic noise. Specifically, 
half-waveplate rotation is limited to modulation frequen- 
cies of tens of Hertz (Hz) at most and boresight rotation 
is limited to much less still. Parallactic angle rotation 
is also a slow signal modulation and is also location de- 
pendent; at South Pole observatories, such as BICEP, 
parallactic angle modulation caused by sky rotation is 
simply not present. 

In this paper we describe the Faraday Rotation Modu- 
lator, a solid-state polarization switch that employs the 
Faraday effect to rapidly modulate polarized light at 
cryogenic temperatures. The FRM is compact, works 
over a large frequency range, and is capable of modula- 
tion rates up to 10 kHz. Specifically, we describe the lab- 
oratory testing and Galactic observations of the FRMs in 
the BICEP experiment during the 2006 observing season. 

2. FRM DESIGN FOR USE IN BICEP 

Figure [T] shows a cross-sectional schematic of a FRM. 
Along the path of incident radiation are two alumina 
cones attached to either side of a ferrite cylinder, pro- 
viding an impedance match to minimize reflection loss. 
The ferrite and cones form a "toothpick" assembly. The 
input and output waveguides preserve polarization ori- 
entation. The ferrite itself is surrounded by a supercon- 
ducting solenoid and the entire toothpick assembly acts 
as a dielectric waveguide, allowing hybrid electric modes 
to propagate. Details of the des i gn an d construction of 
FRMs can be found in lKeatingI (j2009f) . 

Polarization rotation takes place only within the cylin- 
drical ferrite. Bias currents driven through the solenoid 
generate a longitudinal magnetic field within the FRM, 
rotating the polarized input by an angle 



= VIB^ 



(1) 



where V is the Verdet constant, a parameter describing 
the intrinsic polarization properties of the ferrite mate- 
rial, I is the length of the ferrite cylinder, and B is the 
strength of the applied magnetic field. Equation 1 is 
known as the second-order magneto-optical Faraday ef- 
fect. 

The FRMs were optimized for use in BICEP, a ground- 
based CMB polarimeter designed to measure the B-mode 
polarization of the CMB at degree angular scales. BI- 
CEP is an on-axis refracting telescope with a 0.25 m 
aperture capable of scanning in azimuth and elevation 
and rotating about the boresigh t. Full details of th e BI- 
CEP instrument can be found in lYoon et all (|2006[) and 
iTakahashi et al. 1 ()2009f ). Here we provide only a brief 
summary. 

BICEP 's focal plane consists of 49 pairs of orthogo- 
nally oriented pola rization sensitive bolometers (PSBs) 
(| Jones et al. I I2003D cooled to 0.25 K. During the 2006 
observing season, 25 PSB pairs were tuned for band cen- 
ter of 100 GHz (beam size 0.93 degrees) and the other 
24 for 150 GHz (beam size 0.60 degrees). Each pixel is 
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Figure 1. Cross-Sectional sciiematic of a FRM. Polarized radia- 
tion traverses the FRM from riglit to left. The alumina cones serve 
as an impedance match and are attached to the ferrite in the cen- 
ter. The ferrite is supported by dielectric washers (shown in pink). 
A superconducting solenoid surrounds the ferrite and a corrugated 
waveguide surrounds the cones. 



individually coupled to a stack of three corrugated mi- 
crowave feed horns cooled to 4K. The upper edge of the 
frequency pass band is defined by a series of metal mesh 
resonant filters placed in front of the feed horn. The 
lower edge is defined by the waveguide cut-off imposed 
by the horn itself. 

The FRMs were present in the optical path of six PSB 
pairs around the perimeter of the focal plane: three at 
each of 100 GHz and 150 GHz. Figure [2 shows the lo- 
cation of the FRMs in the focal plane of BICEP along 
with a cross section of the FRM positioning in the opti- 
cal path. For the remainder of the paper the FRMs at 
100 GHz will be referred to as lOOA, lOOB, and lOOC 
and at 150 GHz as 150A, 150B, and 150C as labeled in 
Figure a). 

3. INSTRUMENT CHARACTERIZATION 

The FRMs were subjected to a rigorous series of lab 
and field tests to characterize their behavior and choose 
optimal operating parameters for BICEP observations. 
The following section describes the FRM laboratory test- 
ing and results. 

3.1. Rotation Angle and Bias Signal 

FRMs modulate mm-wave signals by rotating the di- 
rection of polarized radiation relative to the axis of po- 
larization of the corresponding bolometers. To mitigate 
systematic offsets, this relative rotation must be well- 
calibrated and time- invariant. Measurements of the ro- 
tation angle as a function of bias current were made by 
placing an aperture-filling wire grid polarizer (resulting 
in a 100 % linearly polarized signal) at the telescope win- 
dow, biasing the FRM with a triangle wave signal, and 
measuring the voltage response of the PSB pair beneath 
the rotator. An example output of this technique is given 
in Figure 131 

As the FRM bias current is increased and decreased, 
the bolometer voltage responds with a sinusoidal wave 
form that "flattens" at the turnarounds when the fer- 
rite begins to saturate magnetically. From this data, the 
relative rotation angle of the FRM, 9, can be calculated 
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Figure 2. (a) The BICEP focal plane for 2006. Each circle repre- 
sents a pixel with a pair of orthogonal bolometers with polarization 
sensitivity axes depicted by -|-. The larger blue circles represent 100 
GHz pairs and the smaller red circles represents 150 GHz pairs. 
The dashed line represents the boresight reference angle of zero. 
The circles that are shaded in gray in boxes around the perimeter 
feature FRMs. (b) The relative location of FRMs within the opti- 
cal path. The FRMs are housed at 4 K between the feedhorn and 
bolometer enclosure. The bolometer signal is demodulated using 
lock-in amplification with phase referenced to the driving current. 



by 



V - (Knax + Knin)/2 



(2) 



where V is the bolometer response voltage and Vi„ax and 
Knin are the maximum and minimum of that voltage, re- 
spectively. Figure [3] shows the voltage response curve as 
a function of 9. FRMs in BICEP can continuously ro- 
tate a polarized signal through a range of approximately 
±80°. Maximum rotation corresponds to a bias current 
of approximately 300 mA. 

The curves shown in red in Figure |3] are for increasing 
bias whereas the curves shown in blue are for decreasing 
bias. The two curves do not overlap due to magnetic 
hysteresis, where a change in magnetization will lag with 
respect to the externally applied field due to configu- 
rational internal forces. Magnetic hysteresis prevents as- 
signing a one-to-one correspondence between current and 
rotation angle, creating a degeneracy in the bias angle for 
values of the bias current where hysteresis is present. Ad- 
ditionally, the degree of hysteresis was shown also to be 
a function of bias frequency, further increasing the com- 
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Figure 3. Raw detector response of a PSB pair (upper two plots) 
and calculated rotation angle of the FRM (lower to plots) is shown 
for a single FRM calibration test. The red and blue curves corre- 
spond to increasing and decreasing FRM bias current respectively. 
FRM hysteresis can clearly be seen in both plots. The hysteresis 
curve separation at zero bias current was between 10° and 25°. 

plexity of demodulation. To avoid the complications of 
magnetic hysteresis the FRMs were biased with a square 
wave signal. This limited the FRM rotation angle to only 
two states, which could be uniquely determined by the 
magnitude and sign of the bias current. 

The amplitude of the bias signal was chosen such that 
the FRMs would have angular modulation of 45 degrees, 
corresponding to a bias current of approximately ±125 
mA. In the reference frame of the bolometer, this is 
equivalent to switching between Stokes +Q and -Q, ef- 
fectively transforming a single PSB into two orthogonal 
detectors. The square wave bias was tuned for a fre- 
quency of 1 Hz, slightly higher than the 1// knee for the 
BICEP t emperature data, gener ally 0.5-1 Hz for a sin- 
gle pixel (jTakahashi et al.l 120091 ) . During unmodulated 
CMB and Galactic half-scans, BICEP shifts the 1// knee 
by scanning the telescope at a scan speed of 2.8° /s. Al- 
though the FRMs are capable of switching at frequencies 
orders of magnitude higher, the FRMs in BICEP were 
limited by the bolometer time constants, which had a 
median of 21 ms across the focal plane. A 1 Hz bias was 
found to be the best compromise between maximizing in- 
tegration time per cycle, which lowers the frequency but 
decreases the noise in demodulation (Section 4.2.3), and 
minimizing the 1// atmospheric fluctuations. 

3.2. Instrumental Polarization 

Instrumental polarization (IP) is spurious polarization 
detected when only unpolarized light is observed. IP can 
result from use of the FRMs if FRM toothpick assembly 
is tilted relative to the optical axis of light to the detector. 
The FRM's IP were measured by placing unpolarized 
aperture-filling 300 K and 77 K sources at the cryostat 
window and biasing the FRMs through their full range 
of rotation angles. The resultant fractional IP is defined 
as 

1 / Vac{300K) \ / 300JC-77X \ 

" 2 V Vdc {300K) - Vdc {77 K) ) V 300K ) 

(3) 
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Figure 4. Dielectric sheet calibration for FRM150A. The top 
panel shows the 1 Hz FRM modulated signal. The bottom panel 
shows the split demodulated timcstreams at ±45 degrees. 

where Vnci'^OOK) is the mean bolometer voltage at 300 
K, VDci77K) is the mean at 77 K, and Vac{SOOK) is 
the peak-to-peak bolometer voltage at 300 K. 

Individual FRM pixels were found to have a repeat- 
able and time-invariant IP signal, with a standard devi- 
ation of approximately six percent (and limited to ten 
percent) over the testing period of several months. The 
scatter of the mean IP between FRM assemblies, how- 
ever, varied considerably. The average value of all FRM 
assemblies was found to be 0.59% with a standard devia- 
tion of 0.43%. The maximum value among all the FRM 
pixels was found to be less than 1.2%, for FRM150A. 

3.3. Rotation Angle Calibration 

Two rotation angle calibration runs were performed 
with the BICEP FRMs. A dielectric sheet calibra- 
tor (DSC) was used based on the design from POLAR 
fO'D eh et al~|[200l iTakahashi et al. 11200811 . The DSC 
consisted of an eccosorb-lined cylinder with a polypropy- 
lene sheet mounted at 45° that served as a beam splitter, 
creating a small polarized signal from input ambient tem- 
perature and sky loads. The direction of the polarized 
signal was determined by the angle of the DSC relative 
to the detector. The DSC was placed above BICEP's 
optical window, and as the boresight was rotated each 
PSB exhibited a sinusoid-like response that varied as a 
function of radial position in the focal plane. Model fits 
using the known properties of the DSC yielded detector 
polarization angles for each PSB in the focal plane. 

Figure 2] shows the tiniestream response of a BICEP 
PSB to DSC calibration with FRM modulation. To ex- 
tract polarization angle, the FRM output was demodu- 
lated into two separate timestreams as described in Sec- 
tion 4.2.3. The result is shown in the bottom panel of 
Figure 4. Each timestream was fit to the DSC model for 
detector angle. Given a nominal PSB orientation angle 
^, derived from DSC tests in the absence of modula- 
tion, the two demodulated timestreams with proper bias 
should yield detector angles V' ± 45°. 

DSC calibration showed proper bias rotation for four 
of the six FRMs. One of the rotator biases was inca- 
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Table 1 

DSC calibration results for all six PSB pairs. The table lists 
the total rotation modulation of the polarization angle. For 
perfect ±45'*'^^ polarization modulation, 90° is expected. 
Polarization modulation in the plus and minus direction is 

one-half of the listed value. The standard deviation is 
calculated across each individual run and averaged over the 
two different rotation angle calibrations. 

FRM PSB Avg Rotation Range ± Standard Deviation (°) 



lOOA 


85.5 ± 7.3 


lOOB 


41.4 ± 12.5 


lOOC 


79.7 ± 7.4 


150A 


90.2 ± 0.1 


150B 


88.2 ± 4.5 


150C 


73.8 ± 18.8 



pable of supplying the currents necessary for full ±45° 
rotation (FRM lOOB) and another was deemed unreli- 
able based on the high standard deviation in its rotation 
angle (FRM 150C). In addition, the PSB pair associated 
with FRM lOOC had excess noise in both the modulated 
and demodulated timestreams, rendering it ineffective for 
celestial polarization studies. 

The results of the DSC calibration are summarized in 
Table [T] for all six PSB pairs. The standard deviations 
listed in the table reflect a true deviation of the rota- 
tion range during bias periods and are not limited by 
measurement readout accuracy. 

4. GALACTIC OBSERVATIONS 
4.1. Observing Strategy 

In 2006, five observations of the Galactic plane were 
made with the FRMs biased. The target fields are shown 
in Figure 5. The first of these observations, referred to 
as the "shallow" observation, is defined by the celestial 
coordinate system (a, 6) = (Right Ascension, Declina- 
tion) as 180° <a < 290° and -70° < 5 < -45° (Figure 
5). The shallow FRM observation was used mainly for 
diagnostic purposes; this was the first test of FRM mod- 
ulation on a celestial source. 

All deep observations focused on a small portion of the 
Galactic plane 238° <a < 248° and -53° < S < -50°. 
The small area was necessary to maximize integration 
time, as the goal of the deep observations was a de- 
tection of polarization from a celestial source. Due to 
their respective positions in the BICEP focal plane, no 
two FRMs could scan this same region of sky simultane- 
ously, so integration time was split evenly between two 
FRMs: FRMIOOA and FRM150A. FRMIOOA was cho- 
sen because it had the highest optical efficiency, lowest 
IP, and a reliable current bias signal. FRM150A was cho- 
sen for its extremely consistent bias signal during DSC 
calibration and its close proximity to FRMIOOA in the fo- 
cal plane. An integration time of approximately 72 hours 
on this deep patch was achieved for each FRM used for 
analysis. 

For all FRM observations, azimuth-elevation raster 
scans were used. At each step in elevation, the telescope 
performed a number of unidirectional scans in azimuth 
(known as "half-scans" ) back and forth across the target 
patch with a scan speed 0.25°/s. The elevation step size 
and number of half-scans at each step were altered be- 
tween the shallow and deep observations. These values 
are summarized in Table [2] along with the other pertinent 
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Figur e 5. 100 GHz FDS M odel 8 Galactic dust emission predic- 
tion ( Finkbeiner et al. Ill999t ) plotted in celestial coordinates on a 
logarithmic scale. The outlined regions show the areas of integra- 
tion for the shallow and deep FRM observations. 

features of each observation. During the deep observa- 
tions, a single boresight rotation angle corresponds to a 
single FRM pixel focused on the Galactic region of inter- 
est. 

4.2. Timestream Processing 

When characterized by Stokes T, Q, and U parame- 
ters, the time domain response of an unmodulated BI- 
CEP PSB as a function of frequency v and direction 51 
is modeled by 

d{t) =Kt® \n{t) +g J dvA,F{v) j dnP{n) 

T+i^(g(cos(2V) + C/sin(2^))]}, (4) 

where V' is the polarization orientation angle of the PSB, 
the parameter e quantifies the amount of cross polar leak- 
age of the PSB, P{fl) is the beam response function, F{i') 
is the spectral response, is the effective antenna area, 
g is the gain factor that converts voltage to temperature, 
n{t) is a term describing total atmospheric and detector 
noise, and Kt is the time domain impulse response of the 
detector. A complete description of the measurements 
and experi mental procedure s used to obtain these terms 
IS given m iTakahashi eraTl (|2009[ ). 

4.2.1. Preliminary Processing 

Data processing begins by deconvolving the complete, 
measured transfer functions for all FRM timestreams. 
The beginning and end of each half-scan is identified 
and periods of non-zero telescope mount acceleration are 
trimmed from the end points. Horizon and celestial bore- 
sight coordinates are calculated by applying a pointing 
model to the raw boresight pointing using the focal plane 
coordinates of the pixel of interest. 

4.2.2. Relative Gains 

Relative detector gains are applied by fitting to the ele- 
vation nods, a period at the beginning of each fixed eleva- 
tion scan when the telescope performs a rounded triangle 
wave motion with an elevation amplitude of 1.2 degrees 
and a duration of approximately 27 seconds. During this 
time, for consistency, the FRMs are modulating the sky 
signal. The responsivity factor for each PSB is found by 
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Table 2 

Details of each FRM Observation, ordered sequentially. Each observation was segmented into 9-hr sections described by the No. 
Sections column. For each of these sections, the rotation angle of the boresight is described. The next two columns describe the 
range of elevation and step size at the boresight center. Azimuth parameters include the length of all azimuthal half-scans, the scan 
speed, and the number of uni-direcational half-scans performed at each step in elevation. 



Date(s) No. Sections Boresight Angle (°) Elevation (°) El Step Az Range (°) Scan Speed (°/s) Half-Scans 



April 21, 2006 5 {315, 315, 180, 180, 0} 55.5-60.5 0.25 ^77 0.25 4 

Oct. 14, 2006 4 {154, 164, -152, -142} 55.5-61.2 0.1 ^12 0.25 22 

Oct. 18, 2006 4 {154, 164, -152, -142} 55.5-61.2 0.1 Sil2 0.25 22 

Oct. 28, 2006 4 {120, -150, -150, 120} 55.5-61.2 0.1 ft!l2 0.25 22 

Oct. 30, 2006 4 {120, -150, -150, 120} 55.5-61.2 0.1 ^12 0.25 22 
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Figure 6. Illustration of time domain demodulation. The top 
plot shows the bias current as a function of time. Points in pink 
are positively biased, corresponding to a -(-45° rotation. Points in 
blue are negatively biased, corresponding to —45° rotation. The 
black crosses are the midpoints of each bias plateau. The middle 
plot shows the corresponding points mapped onto the bolometer 
timestream. The red crosses are points of gradual switching re- 
sultant from bolometer time constants, and are cut from analysis. 
The remaining points are averaged together for each half-plateau 
and are stored as two new demodulated timestreams as shown in 
the lowest panel. The pink filled circles comprise the demodu- 
lated +45° timestream, whereas the filled blue circles correspond 
to —45° timestreams. 

fitting the detector voltage versus the changing line of at- 
mospheric sight, given by the cosecant of elevation mul- 
tiplied by the scale height. Before applying the relative 
gains, each PSB gain factor is weighted by the average 
of the PSB pair during the scan set. 

4.2.3. Demodulation 

Although many demodulation techniques were consid- 
ered, time domain demodulation was used due to its 
simple and robust properties and accurate polarization 
reconstruction. Figure 6 is a graphical example of the 
timestream demodulation using FRM 150A. 

The first step of demodulation was to clean the bias 
current signals of software spikes and glitches via nearest- 
neighbor interpolation. The current signal then dis- 
played fluctuations on the < 1% level, and so perfect 
±45° rotation was assumed for positive and negative val- 
ues of the bias, respectively. 

The current biases were then used as the lock-in phase 
reference signals for demodulation of the corresponding 
bolometer timestreams. All points at which the cur- 
rent bias transitioned across zero were identified and 



the median number of points between transitions was 
computed. Any region where the switching period dif- 
fered from the median by more than 10%, along with 
one plateau region on either side, was identified and the 
corresponding region in the bolometer timestream was 
excluded from analysis. This resulted in a loss of < 1% 
of the usable timestream data as the bias signal exhibited 
very few deviations from normal periodicity. 

The midway point between each two adjacent tran- 
sitions in the bias signal was identified, dividing each 
"bias plateau" into two halves. The cleaned bias signal 
and its segmentation are shown in the top panel of Fig- 
ure 6. The average value of the corresponding bolometer 
timestream, along with other pertinent variables (point- 
ing, time, etc.), for each half was computed, less the two 
to four points closest to the transition to account for the 
non-instantaneous bolometer response. Although decon- 
volution using the measured time constant does recover 
a small number of points along each transition - with- 
out deconvolution, a minimum of six points were cut on 
each side ~ the bolometer response to bias switching is 
still imperfect and a fraction of the data surrounding 
each transition must still be removed from analysis. The 
number of points cut from each side of the transition was 
derived independently for each bolometer using iterative 
fits to the elevation nods. To maximize the signal-to- 
noise, the minimum number of data points was cut that 
yielded a convergent fit to the elevation nod signal. 

The bolometer timestream averaged values were then 
separated based on the sign of the corresponding bias sig- 
nal. Essentially, the timestream of a single bolometer at 
polarization angle V're/ has been demodulated and down- 
sampled into two demodulated bolometer timestreams at 
V're/ + 45° and V're/ — 45°. This is shown in the last panel 
of Figure 6. The sum and difference of the two demod- 
ulated timestreams from each individual PSB was then 
calculated. 

4.2.4. Masking and Filtering 

A DC-ofFset was removed from both the sum and differ- 
ence timestreams via mean-subtraction of each azimuthal 
half-scan. An obvious distortion arises when fits for the 
mean include the strong Galactic signal, so the Galaxy 
was masked from the fits. The top two panels of Fig- 
ure [7] show an example of such distortion and the lack of 
distortion in the masked fit. 

BICEP CMB ha lf-scans are filtered using a third 
order polynomial ()Chiang et al. I I2O1O0 and BICEP 
Galactic half-scans u tilize a second order polynomial 
(|Bierman et al. Il2011f) in order to accurately remove at- 
mospheric 1// noise. For the deep FRM observations, 
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the length of the half-scans across the galaxy prohib- 
ited such a filtering scheme. A mask large enough to 
remove all Galactic signal from the fits resulted in the 
majority of half-scans beginning or terminating on the 
masked region. Interpolated fits that are not constrained 
by data on both sides of the mask tend to diverge and 
no longer reasonably approximate noise within the mask. 
Conversely, smaller masks avoid this problem but leave 
residual Galactic signal for the polynomial fit, causing 
higher order polynomial fits to remove true Galactic sig- 
nal instead of only noise. These results are shown in the 
last two panels of Figure [T] While mean subtraction does 
not remove large-scale noise from the data, FRM modu- 
lation mitigates some of the 1 // atmospheric fiuctuations 
that would otherwise need to be filtered. 

The optimal mask cuts maximal Galactic signal from 
the fits while retaining sufficient off-Galactic data to 
yield statistically significant fits for the mean. Given the 
Galactic coordinates defined by (&, I) = (latitude, longi- 
tude), the radius of the masked region was reduced from 
I 6 |< 3° (the maximal possible mask for the deep inte- 
gration scans) until the DC offset fits converged for the 
majority of half-scans. The optimal mask was found to 
be I 6 |< 1.5°. 

Following DC offset subtraction, residual high- 
frequency power from demodulation was removed via a 
Butterworth low pass filter at 0.6 Hz. A data quality 
check was then performed where any half-scan with a sig- 
nal spike more than 7 times the standard deviation was 
removed from analysis. These regions comprised < 1% 
of the total usable data. 

4.2.5. Instrumental Polarization Removal 

Temperature to polarization leakage from a bright 
source such as the Galaxy can cause significant artificial 
structure in polarization. Laboratory testing revealed 
that FRM150A showed approximately 1% IP, the highest 
of any of the devices; a value that could distort Galactic 
polarization maps. To remove this temperature to polar- 
ization leakage, the IP value for each half-scan was cal- 
culated using a fit to the preceding elevation nod. The 
derived fractional IP value was then multiplied by the 
half-scan sum data and subtracted from the difference 
data. 

Figure [8] shows an example of IP leakage derived from 
elevation nods. Because the atmosphere is unpolar- 
ized, the elevation nods should cause a change in the 
summed (but not differenced) data for each demodu- 
lated PSB. The data from each elevation nod for each 
PSB was demodulated, summed/differenced, mean sub- 
tracted, and low-pass filtered. The summed and differ- 
enced timestreams are both fit to the cosecant of eleva- 
tion to derive average responses to the changing line of 
sight. Average IP values were computed by taking the 
mean of the absolute value of the difference divided by 
the sum across all scans and all runs. The average values 
at both frequencies across all observations was found to 
be IP = 0.21% and IP = 0.88% at 100 GHz and 150 
GHz respectively. 

To confirm that the elevation nods revealed IP and not 
some other source of spurious polarization, the fit values 
for both the sum and difference were plotted as a function 
of elevation. The summed data increased with elevation 
due to greater atmospheric loading on the PSBs, and the 
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Figure 7. Summed detector timestreams for FRM deep integra- 
tion Galactic half-scans witli different filtering methods. The top 
plot shows the timestream (blue), DC offset fit (red), and the resul- 
tant filtered timestream (black). Note the severe distortion from 
fitting the Galactic region. The second plot shows the same data 
with Galactic intensity removed from the fit. The boundaries of 
the mask are shown by black dotted vertical lines. The DC offset 
is subtracted from the masked region, leading to minimal artifacts. 
The third plot shows data that terminates on the masked Galactic 
region, fit with a higher-order polynomial; the interpolation di- 
verges through the Galactic plane. The final plot illustrates the 
effects of subtracting a higher order polynomial fit from data with 
an insufficient mask to completely remove Galactic signal. 

differenced data was shown to scale proportionally to the 
sum. The correlation between summed and differenced 
data indicate true IP. 

4.2.6. Map Making 

After timestream data processing was complete. Stokes 
T, Q, and U were derived following the formalism in 
iJones et al. I ^007). Temperature data were obtained 
by binning the filtered detector timestreams into map 
pixels, p, as 



E 



where the indices i and j denote the PSB channel num- 
ber and timestream sample number respectively. The 
variable is the sum demodulated timestream, n is the 
number of PSBs, and w+ is the weight assigned to each 
demodulated sum. Stokes Q and U were calculated by 
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Figure 8. From top to bottom: Elevation angle, FRM150A de- 
modulated sum, and FRM150A demodulated difference data for 
an elevation nod versus time for FRM 150A. The mean has been 
subtracted from both the sum and the difference data so that small 
changes can be seen. The difference data is approximately 1%, of 
the sum data, indicative of instrumental polarization. 
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where w~ is the weight assigned to the demodulated 
timestream difference and a and (3 are the differenced 
orientation factors 
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The factors of ±45 degrees come from the Faraday ro- 
tation of the nominal polarization angle of the i*'* PSB. 
The same correction factor describing cross polar leakage, 



^, is applied to both terms and Pu. Equa- 



tion 5 is degenerate for a single FRM PSB at time j. To 
break the degeneracy, the same sky pixel p is observed 
by the same modulated detector at a different polariza- 
tion angle V'lj • With more than one observing angle the 
off-diagonal elements of Equation 5, aijPij, average to 
zero and the matrix can be inverted to solve for Q and 
U. 

For each half-scan the inverse of the variance of the pair 
sum and difference timestreams is used as the weight . 
Total integration time is also computed for each pixel. 

4.2.7. Absolute Calibration 

To relate detector units to CMB units, BICEP mea- 
surements are cross-correlated with WMAP temperature 
maps to derive an abso lute gain calibration . A complete 
description is given in (IChiang et al. l[2010l) . 

This same comparison cannot be made for the FRMs, 
as the FRMs were never biased and observing during 
CMB observations. Instead, FRM detectors were cali- 
brated via comparison with the calibrated 3-year BICEP 
temperature maps of the bright arm of the Galaxy. The 
absolute calibration factor was computed as the slope 
of a linear fit to the pixel-pixel scatter plot of the two 
maps. The variance of the difference between the fit 
data and the BICEP data is used to compute the un- 
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Figure 9. Integration time and temperature maps from the shal- 
low FRM observation on April 21, 2006. The temperature of the 
Galactic plane can clearly be seen above the residual noise. The 
large scale noise features at right ascension 180° — 200° and decli- 
nation of ft! —63° are an artifact of masking and filtering (Section 
4.2.4). Although the non-FRM PSBs are observing during this 
shallow observation, the slow scan speed places their observing 
band well below the 1// knee, leading to excessive noise. For this 
reason, these maps are intentionally omitted from further analysis. 

certainty. This method yields values of 1.1288 ± 0.0023 
and 1.0011 ± 0.0017 times the nominal absolute calibra- 
tion factors for BICEP at 100 and 150 GHz respectively. 

5. MAP RESULTS 

Figure [9] shows both the integration time and Galac- 
tic temperature maps derived from the shallow FRM 
observation at 100 GHz and 150 GHz. Because the 
shallow observation featured long sweeps of the Galaxy, 
half-scans were filtered using a second order polynomial 
rather than mean subtraction. T he maps were binned 
with Healpix (jCorski et al. Il2005[ ) pixelization 0.25°. A 
a = 0.5° Gaussian smoothing function has been applied 
to all maps. 

Due to glitches in the rotator bias signals, only the 
first and last 9 hour sections resulted in usable data, 
a total of approximately seventeen hours of integration 
time including data cuts. The deepest integration time 
was « 15s per pixel, with an average of 4.1 and 5.4 s 
per pixel at 100 and 150 GHz, respectively. Data from 
5 out of 6 of the FRM pixels (two at 100 and three at 
150) were used for these maps; one was omitted due to 
biasing problems. 

The temperature maps from this observation reveal 
bright Galactic emission confined mainly to | 6 |< 1°, 
with both large scale features and compact sources. The 
intensity is greater at 150 GHz than at 100 GHz, as 
predicted by models of dust in the interstellar medium 
(|Finkbeiner et al.1ll999[ ), and consistent with other re- 
cent mm- wave exper iments (jCulverhouse et alTI 120101 : 
iBierman et al. Il2011[ ). Although the off-Galactic por- 
tions of the map show relatively large noise fluctuations 
due to the low integration time, the Galactic signal domi- 
nates the noise for both bands. This is, however, not true 
for the much fainter polarization signal (not pictured); a 
detection of statistically-significant polarization was not 
seen in this observation. The large scale noise features 
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that differ between the 100 GHz and 150 GHz maps re- 
sult from a combination of 1// atmospheric noise and 
beam smoothing. 

The left two cohimns of Figure 10 show integration 
time as well as the Galactic T, Q, U, polarization vec- 
tor, and Q and U difference maps at both 100 GHz 
and 150 GHz from the FRM deep observations. Only 
two FRMs, lOOA and 150A, were used to accumulate 
the m aps. Maps are binned into Healpix (jGorski et al. I 
[2005h 0.1° pixels and high -pass filtered via DC offset sub- 
traction (section 4.1.1). K a — 0.2° Gaussian smoothing 
kernel has been applied to the maps at both frequencies. 

The integration time in this region averages 77.7 and 
67.2 seconds per map pixel at 100 and 150 GHz respec- 
tively. The T, Q, U maps all show an increase in signal 
from 100 to 150 GHz. At both frequencies, large-scale 
positive Q and negative U polarization can be seen along 
the plane of the Galaxy. These polarization features are 
generally confined to within | & |< 1.5° of the center of 
the Galactic plane. The polarized portions of the maps 
are coincident with the largest magnitude intensity in 
Declination and slightly offset towards lower Right As- 
cension. At both frequencies, the magnitude of the signal 
in Stokes U is approximately five times larger than the 
signal in Stokes Q. The average signal in Stokes U con- 
fined to I 6 |< 1.5° is given by -16.2 fiK and -48.4 jiK at 
100 and 150 GHz respectively. Given the error bars de- 
fined by difference maps (Section 5.2), this corresponds 
to an average S/N ratio of 1 and 3 in each smoothed map 
pixel for the average galactic polarization signal seen at 
100 and 150 GHz respectively. 

BICEP observations of the same region were made dur- 
ing the 2006 - 2008 BICEP observing seasons and com- 
pared to the FRM data. The BICEP Galactic scan strat- 
egy uses long sweeping scans across the entire Galactic 
plane at a step size of 0.25°, prohibiting a direct com- 
parison of a single non-FRM PSB pair to a single FRM 
PSB pair; for a single non-FRM PSB pair, the integra- 
tion time is too low and the number of boresight an- 
gles is insufficient to accumulate polarization. As such, 
maps were accumulated using all non-fiagged PSB pairs 
in the focal plane. PSB pairs containing FRMs (among 
others) were flagged and excluded from analysis. The 
data were taken over all three BICEP observing seasons, 
totaling 763 observing hours for the entire shallow ob- 
servation region. Non-FRM data were processed using 
the same analysis pipeline as the FRM data with a few 
exceptions. The initial deconvolution was applied with 
a low-pass filter at 5 Hz. Following deconvolution, the 
data were downsampled to 10 Hz before preliminary pro- 
cessing and relative gain calibration. The demodulation 
step was omitted; the sum and difference data were taken 
between the two PSBs within each pair. The data were 
then subjected to the same masking and filtering strat- 
egy as the FRM pixels, with the additional complication 
that all data outside the region of deep integration were 
also masked from the fits. This is to keep large por- 
tions of off-Galactic data from weighting the DC offset 
subtraction, which was not possible with the FRM deep 
integration scan strategy. Maps were accumulated using 
the same formalism as Section 4.2.6, except Equations 
(6) and (7) are modified to 



Pi] = 7i(sin[2('0i+ij)] - sin[2(V'y)] 



(9) 



where i and i + 1 are the indices of the two PSBs within 
the pair that have been subtracted. The binning and 
smoothing of the maps is identical to the FRM deep ob- 
servations. 

The right two columns of Figure 10 show the integra- 
tion time as well as Galactic T, Q, U, polarization vec- 
tor, and Q and U difference maps at both 100 GHz and 
150 GHz from the BICEP observations accumulated us- 
ing the non-flagged focal plane PSBs. A total integration 
time of approximately 324 hours and 194 hours was taken 
at 100 and 150 GHz respectively, with a mean integration 
time of 412 seconds per pixel at 100 GHz and 246 seconds 
per pixel at 150 GHz. The maps exhibit very good agree- 
ment with the FRM maps in both temperature and po- 
larization at both frequencies. As evidenced by the polar- 
ization vector maps in Figure [TOl both data sets show the 
same strong Galactic polarization, roughly perpendicular 
to the Galactic plane, with the highest polarized signal 
in the upper left quadrant. Both the FRM and non- 
FRM BICEP maps exhibit a diminution of polarization 
across the Galaxy with increasing RA. At each frequency 
and for all polarization maps, the contribution of Stokes 
U dominates the polarized signal. Both the FRM and 
non-FRM BICEP maps are consistent with the polarized 
emission expectations based on Galactic plane maps pub- 
lished previo usly by BICEP and other recent mm-wave 
experimen ts (iHildebrand et ah lll999l:lMeiia et al 



2005 



2010; 



■ 7j(cos[2(V'j+ij)] - cos[2(V'y)]) 



(8) 



Mivillc D eschenes et al. l l2008HCulverhouse et al. 
iBierman et al. 11201 It ). 

The polarization fraction is found by calculating the 
median value of the quotient of polarization and tem- 
perature across all map pixels where | h |< 1.5°. The 
variance in this same region is used to compute the un- 
certainty on these values. The FRM maps yield a po- 
larization fraction of 1.32% ± 2.17% and 2.36% ± 0.21% 
at 100 and 150 GHz respectively. The results derived 
from BICEP 's non-FRM maps show strong agreement: 
1.53% ± 0.61% at 100 GHz and 2.31% ± 0.02% at 150 
GHz. For both the FRM and non-FRM BICEP maps, 
150 GHz displays little variance in polarization fraction 
across the galactic plane, whereas the 100 GHz variance 
is much higher. 

5.1. Noise and Systematics 

The deep FRM observations were analyzed for noise 
and other transient issues that may have resulted from 
the use of the devices. The results of the FRM observa- 
tions were also subjected to many self-consistency checks 
to in order to verify the accuracy of the data presented 
here. Due to the low signal-to-noise achieved for the 
shallow observations, only the deep observations are sub- 
jected to the difference map and noise analysis. 



5.2. Difference Maps 

Difference maps, in which the map data is split in half 
and differenced, were used to check the self-consistency 
of the Galactic maps. Although some residual signal may 
remain due to timestream filtering effects, the expected 
signal of the maps is nearly zero and all large scale Galac- 
tic temperature and polarization signals should vanish. 
Statistical polarization errors are quantified by taking the 
standard deviation for all Q and U map pixels at both 
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Figure 10. FRM and BICEP maps of the deep integration region. From top to bottom, each column displays a single plot of temperature 
with polarization vectors (T & P), integration time (Tint), Stokes Q, Stokes U, and Q and U scan direction difference maps (Diff Q and 
Diff U). All maps are displayed in /iK except for integration time, which is displayed in seconds. The first two columns are the FRM deep 
observation results at 100 and 150 GHz. The last two columns show maps made with the nominal BICEP observations using non-FRM 
PSBs at 100 and 150 GHz. 
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Figure 11. T, Q, and U differenced data from individual detec- 
tors within a single PSB pair. The Q and U scales are identical 
to Figure 7 for comparison, whereas the temperature scale is re- 
duced to ±500 fiK. Large scale features in both polarization and 
temperature are no longer apparent in the maps. 

100 and 150 GHz. These results are summarized in Ta- 
ble [3] For consistency with the map results, a. a = 0.2° 
Gaussian smoothing function is applied to all maps. As 
such, the rms map noise for each pixel in this Table [3] is 
quoted for an effective pixel size of 0.24 deg^. 

Three separate divisions of the data were analyzed: 
scan-direction, PSB pair, and observation. The scan- 
direction split, where the two data sets are separated 
based on half-scan direction, can generate noise based 
on thermal instability at half-scan endpoints. The scan- 
direction difference maps are shown in the last two rows 
of Figure (TU] PSB pair difference maps subtract the 
maps accumulated from individual detectors within a 
pair. This difference map is unique to FRM analysis, 
as a FRM turns a single BICEP PSB into a polarime- 
ter, allowing T, Q, and U maps to be accumulated for 
individual detectors within a PSB pair. This difference 
map is perhaps the most robust data quality test for the 
FRMs as it is sensitive to many factors including rela- 
tive gain mismatches, demodulation errors between de- 
tectors, and thermal stability. The PSB pair difference 
maps are shown in Figure [TTJ Finally, observation dif- 
ference maps subtract maps accumulated from odd and 
even numbered observing runs. This jackknife probes 
for sensitivity to weather changes. The FRM maps show 
the highest sensitivity to changes in weather, especially 
in the 150 GHz PSB. The main source of the elevated 
signal at 150 GHz is a very high noise contribution at 
the lowest declination, where the least integration time 
occurs. 

Difference maps of the two data sets, the FRM and 
non-FRM BICEP, are displayed in Figure [T^l The rms 
map noise per pixel is given by 20.6 /xK rms/pixel at 100 
GHz and 35.6 /zK rms/pixel at 150 GHz. This differ- 
encing should remove all true sky signal and leave only 
uncorrelated noise, which should combine as the quadra- 
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Figure 12. Maps of the difference between the FRM and non- 
FRM BICEP observations. Galactic signal, in both temperature 
and polarization, has been removed. This is particularly evident in 
celestial U, where both the FRM and non-FRM BICEP maps origi- 
nally exhibited significant polarized structure. No similar structure 
is evident in the Stokes U difference map. 

ture sum of the noise in the two individual maps. To 
estimate the map noise in the non-FRM BICEP map, 
the scan direction difference map was utilized, yielding 
12.5 and 12.0 per effective pixel at 100 and 150 GHz 
respectively. Given the map noise values in Table [TT] 
for the FRM map, the difference map noise is consistent 
with the quadrature sum of the map noise of the FRM 
and non-FRM BICEP difference maps. 

Although many properties of the FRMs were charac- 
terized during initial testing and deployment, the lack of 
time for field testing prevented a complete characteriza- 
tion for the FRMs and is beyond the scope of this pa- 
per. Specifically, properties such as cross-polar leakage 
and polarization inefficiency were not studied in-depth, 
and it was assumed throughout this analysis that the 
FRMs behaved as ideal devices in these respects. Fur- 
ther, the DSC calibration revealed an imperfect bias sig- 
nal for most of the FRMs, but for the deep observations 
it was assumed that the FRMs accomplished an ideal ro- 
tation of ±45° to the nominal angle of the coupled PSBs. 
Deviations from this perfect modulation would result in 
miscalibrated detector angles and reduced polarization 
efficiency, which would alter the polarization signal in 
the Q and U maps. 

Although the FRM maps are absolutely calibrated 
via non-FRM BICEP maps, the calibration is done us- 
ing temperature alone. Deviations from the (assumed 
perfect) polarization properties of the FRMs (described 
above) would result in marked differences between the 
FRM and non-FRM Q and U difference maps. The fact 
that these maps are statistically consistent with each 
other limits the presence of FRM non-idealities to the 
percent level. 

5.2.1. Map Noise 
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Figure 13. Plot of rms map noise for the scan-direction difference 
maps as a function of the portion of total integration time used. 
At each frequency, two lines are displayed representing the rms 
difference map noise from individual PSB and the FSB pair. The 
thinner dashed lines are a guide to the eye exhibiting square-root 
of integration time dependence; they are plotted in the same color 
as their corresponding solid lines. At 100 GHz, the square-root of 
integration time scaling is obvious, whereas for a single PSB at 150 
GHz there is no such dependence. In addition, the 100 GHz PSB 
pair scales as the square-root of the number of detectors (\/2) from 
the individual PSBs to the pair. At 150 GHz, the \/2 reduction in 
noise is not seen. 

Map noise in FRM observations 2-5 was quantified via 
the scan-direction difference map. If the map was not 
sensitive to systematic effects (e.g. - systematics are not 
present), then the noise should be Gaussian white and 
integrate down with the square-root of the integration 
time and the square-root of the number of detectors. To 
vahdate this assumption for these maps, the data were 
reprocessed several times, each time removing a different 
number of points from the transitions during demodula- 
tion. Maps were accumulated using a range of cuts from 
four points (two on either side of the transition) to twelve 
points (six on either side of the transition) . The results 
are shown in Figure 1131 The data exhibit a decrease in 
rms map noise with a \ft dependence based on the data 
cuts. Additionally, at 100 GHz, the data show a \f2 in- 
crease in the map noise when accumulating the difference 
maps with a single PSB as opposed to the pair. At 150 
GHz, however, the \/2 dependence on PSB number is not 
observed. 

Histograms of the noise distribution were also com- 
puted using the unsmoothed difference maps at both fre- 
quencies and are shown in Figure 1141 Only map pixels 
with an integration time of 50 seconds or more are used 
in this analysis to avoid pixels with low integration time 
weighting the result. The histograms show that the am- 
plitude distribution of the noise is roughly Gaussian and 
distributed about zero within one standard deviation of 
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Figure 14. Histograms of the noise at 100 GHz (left) and 150 GHz 
(right) in the unsmoothed scan direction jackknife maps for the 
FRM deep observations. Gaussian fits using the derived mean and 
standard deviation at each frequency are shown for comparison. 

the mean. The mean ± standard deviation are given by 
-4.1 ± 75.3 ^fiC and - 1 1.1± 99.4 /iif at 100 GHz and 150 
GHz respectively. 

Assuming these maps are Gaussian white noise- 
dominated, detector noise was quantified by multiplying 
the Q and U pixel values by the square-root of the inte- 
gration time per pixel and taking the standard deviation 
across the maps. This yielded an average instantaneous 
(i.e., single Stokes parameter) "NEQ per feed" of 420 
[iK^ and 501 ijlK^ for 100 GHz and 150 GHz respec- 
tiv ely. The value at 100 GHz is comparable to th e values 
in (|Takahashi et al. 1120091 iBierman et 3111201 It ), though 
the value at 150 GHz is elevated by a factor of 1.6. The 
source of the excess noise at 150 GHz is unknown, though 
the lack of dependence when doubling the number of 
PSBs indicates that systematic effects, such as detector 
cross-talk, can not be ruled out. 

6. RESULTS FROM UNBIASED FRMS 

For the majority of the 2006 observing season, the 
FRMs remained unbiased in the BICEP focal plane. 
Without modulation, the PSB pairs coupled to the FRMs 
functioned as nominal BICEP PSBs, where the sum and 
difference of the pair were used to calculate temperature 
and polarization respectively. Using 2006 Galactic scans, 
noise was computed for all FRM feeds to analyze whether 
the presence of the FRMs in the optical path introduced 
additional noise within the detectors. 

To compute the noise, the spectral power distributions 
for detector sums and differences for each calibrated half- 
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Table 3 

The rms map noise for each of the three difference maps analyzed for the deep 
FRM observations. The rms is quoted for an effective smoothed pixel size of 
0.24 deg^. For comparison, the Galactic U signal for the FRM maps was 
found to be -16.2 and -48.4 at 100 and 150 GHz respectively. 



Difference Map Type 100 GHz rms (^jK/pixel) 150 GHz rms (^jK/pixel) 

Scan Direction 12.3 20.5 

PSB Pair 13.4 12.8 

Observation 22.6 45.0 
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Frequency (Hz) 

Figure 15. Average power spectral distributions calculated for 
the pair-sum and pair-differenced FRM PSBs at 100 and 150 GHz 
during Galactic observations in 2006. Data for the plots were taken 
while the FRMs were unbiased in the BICEP focal plane. The 
second order polynomial filtering can be seen at 0.05 Hz, where 
the pair-difference noise is lower than the white noise floor. 

scan is calculated and averaged for all FRM pixels at each 
band. This is done after masking, filtering, and sum- 
differencing each PSB pair. FigurefTSlshows the results of 
this analysis. The pair-sum data exhibits increasing 1// 
contamination from 100 to 150 GHz. Pair differencing re- 
moves this contamination, resulting in nearly white noise 
above 0.1 Hz. Averaging the pair-difference periodogram 
between 0.1 and 1.0 Hz gives a Noise Equivalent Temper- 
ature (NET) per detector of 620 and 430 ^iKy^ for FRM 
feed s at 100 and 150 GHz respectively. Previous stud- 
ies (|Takahashi et al. 1 120091 : iBierman et al. Il2011[ ) char- 
acterizing the pixels in BICEP (which do not include the 
FRMs) give BICEP NET values of 525 fJ-KyG at 100 GHz 
and 450 fiKy/s at 150 GHz. The elevated noise at 100 
GHz was found to stem from the PSB pair corresponding 
to a single FRM (lOOC). When removed, the FRM NET 
at 100 GHz is 520 ^Ky/s. Both bands show excellent 
agreement with the quantities derived from observation 
quality pixels, demonstrating that unbiased FRMs do not 
introduce any excess noise into the detectors. 

7. CONCLUSIONS 

In this paper, we report on the first detection of Galac- 
tic polarization using fast, active optical polarization 
modulation. The observed polarized signal is consistent 
with both BICEP 's and other recent experiments' maps 
of Galactic polarization. Although observation and test- 
ing time limited the scope of the comparison that could 
be made between the nominal BICEP instrument and 
the BICEP instrument modulated by FRMs, the FRMs 
have been shown to be effective solid-state polarization 
modulators for both laboratory and celestial polarime- 
try applications. In particular, the 100 GHz FRM was 
found to exhibit equivalent noise to similar BICEP PSB 
pairs without FRMs, but expands the parameter space 
for the instrument's scan strategy, relaxing constraints 
on scan speed and potentially increasing the observed 
sky fraction. Ultimately, the choice of whether to use 
a fast active modulator such as a FRM depends upon 
the details and constraints of the experiment. Devices 
similar to the FRMs are currently being explo red for use 
in several experiments ((Tartari et al. ..2009; .Gault et al. I 



1 201 J) . FRM technology could be applied wherever fast 
solid-state polarization modulation is called for. 

FRMs are a promising technology due to their design 
flexibility, low systematic polarization, large bandwidth, 
and ability to be used over a large frequency range. Fur- 
ther, FRMs have been shown to be functional in both 
modulating and non-modulating modes without adding 
noise, making them a flexible option for a mm-wave po- 
larization modulation. FRMs also have potential appli- 
cation as modulators for cm and sub-mm wavelength po- 
larimeters. 
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